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A general  formulation  for  aperture  problems  is  given  in  terms  of  the  method 
of  moments.  It  applies  to  any  two  regions  isolated  except  for  coupling  through  the 
aperture.  The  aperture  characteristics  are  expressed  in  terms  of  two  aperture 
admittance  matrices,  one  for  each  region.  The  admittance  matrix  for  one  region  is 
Independent  of  the  other  region,  and  hence  can  be  used  for  any  problem  involving 
that  region  and  aperture.  The  solution  can  be  represented  by  two  generalized 
n-port  networks  connected  in  parallel  with  current  sources.  The  current  sources 
are  related  to  the  tangential  magnetic  field  which  exists  over  the  aperture  region 
when  the  aperture  is  closed  by  an  electric  conductor.  Explicit  formulations  are 
given  for  two  problems, that  of  an  aperture  in  a conducting  plane  with  plane-wave 
excitation,  and  that  of  a vTavegulde  feeding  an  .iperture  in  a conducting  plane.  ' 
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I.  INTROWJCTI(» 


The  general  probleo  of  coupling  through  apertures  has  many  specific 
applications,  such  as  apertures  in  a conducting  screen,  waveguide-fed 
apertures,  cavity-fed  apertures,  wavegulde-to-wavegulde  coupling,  wavegulde- 
to-cavity  coupling  and  cavity-to~cavity  coupling.  The  literature  on  these 
problems  is  extensive.  Many  books,  of  which  [1]  to  [3]  are  typical,  dis- 
cuss the  problem  and  give  references  to  some  of  thu  literature. 

This  report  formulates  the  problem  in  terms  of  a moment  solution  of 
the  operator  equation.  An  application  of  the  equivalence  principle  sepa- 
rates the  problem  Into  two  parts,  namely,  the  regions  on  each  side  of  the 
aperture.  The  only  coupling  Is  through  the  aperture,  whose  characteristics 
can  be  expressed  by  aperture  admittance  matrices,  one  for  each  region.  These 
admittance  matrices  depend  only  on  the  region  being  considered,  being  Inde- 
pendent of  the  other  region.  The  aperture  coupling  is  then  expressible  as 
the  sum  of  the  two  Independent  aperture  admittance  matrices,  with  source 
terms  related  to  the  Incident  magnetic  field.  This  result  can  be  Inter- 
preted In  terms  of  generalized  networks  as  two  n-port  networks  connected  in 
parallel  with  current  sources.  The  resultant  solution  Is  equivalent  to  an 
n-term  variational  solution. 

Since  the  problem  Is  divided  Into  two  mutually  exclusive  parts,  one 
can  separately  solve  a few  canonical  problems,  such  as  apertures  In  con- 
ducting screens.  In  waveguides,  and  in  cavities,  and  then  combine  them  In 
the  various  permutations  mentioned  above.  Computer  programs  can  be  de- 
veloped for  treating  broad  classes  of  canonical  problems,  such  as  apertures 
of  arbitrary  shape  in  conducting  planes,  in  square  waveguides,  and  in  rec- 
tang'jiar  cavities.  Such  programs  can  then  serve  as  broad  and  versatile 

[1]  R.  E,  Collin,  "Field  Theory  of  Guided  Waves,"  McGraw-Hill  Bock  Company, 
New  York,  1960. 

[2]  R.  F.  Harrington,  "Tine-Harmonic  Electromagnetic  Fields,"  McGraw-Hill 
Book  Company,  New  York,  1961. 

[3]  N.  Ainltay,  V.  Tallndo,  C.  P.  Wu,  "Theory  and  Analysis  of  Phased  Array 
Antennas,"  Wlley-Interscience,  New  York,  1972. 
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1.  The  general  problem  of  two  regions 
coupled  by  an  aperture. 


tools  for  deslgaing  electromagnetic  networks  with  aperture  coupling. 


II.  GENERAL  FORMULATION 

Figure  1 represents  the  general  problem  of  aperture  coupling  between 
two  regions,  called  region  a and  region  b.  In  region  a there  are  impressed 
sources  j4,  and  region  b is  assumed  source  free.  The  more  general  case 
of  sources  in  both  region  a and  region  b can  be  treated  as  the  superposition 
of  two  problems,  one  with  sources  in  r:^glon  a only,  plus  one  with  sources  in 
region  b only.  Each  region  of  Fig.  1 is  showii  to  be  bounded  by  an  electric 
conductor,  although  other  types  of  electromagnetic  isolation  may  be  used. 
Region  a is  shown  closeri  and  region  b is  shown  open  to  infinity,  although  each 
region  may  be  open  or  closed.  The  equivalence  principle  [2,  Sec.  3-5]  is  used 
to  divide  the  problem  into  two  equivalent  problems,  as  shown  in  Fig.  2.  In 
region  a,  the  field  is  produced  by  the  sources  J^,  M^,  plus  the  equivalent 


magnetic  current 


M - n X E 


over  the  aperture  region,  with  the  aperture  covered  by  an  electric  con- 
ductor. In  region  b,  the  field  is  produced  by  the  equivalent  magnetic  cur- 
rent -M  over  the  aperture  region,  with  the  aperture  covered  by  an  electric 
conductor.  The  fact  that  the  equivalent  current  in  region  a is  -HM  and  that 
in  region  b is  -M  ensures  that  the  tangential  component  of  electric  field 
is  continuous  across  the  aperture.  The  remaining  boundary  condition  to  be 
applied  is  continuity  of  the  tangential  component  of  magnetic  field  across 
the  aperture. 

The  tangential  component  of  magnetic  field  in  region  a over  the  aper- 
ture,  denoted  H.,  is  the  sum  of  that  due  to  the  impressed  sources,  denoted 
M^,  plus  that  due  to  the  equivalent  sources  M,  denoted  that  is 

«) 

^ a 

Note  that  H and  If  (M)  are  both  computed  with  a conductor  covering  the 
aperture.  A similar  equation  holds  for  region  b,  except  that  the  equiva- 
lent sources  -M^ are  the  only  sources.  Hence,  the  tangential  component  of 
magnetic  field  in  region  b over  the  aperture  is 

^ = - h!?(M)  (3) 

(wt  ^ >*»• 
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where  is  coeluted  with  a conductor  covering  the  aperture.  The 

last  equality  In  (3)  Is  a consequence  of  the  linearity  of  the  h|^  operator. 
The  true  solution  Is  obtained  when  of  (2)  equals  of  (3) , or 


H*((p  - -hJ 


(4) 


This  is  the  basic  operator  equation  for  determining  the  equivalent  mag- 
netic current  M. 

If  (4)  were  satisfied  exactly,  we  would  have  the  true  solution. 

We  use  the  method  of  moments  [4]  to  obtain  an  approximate  solution.  De- 
fine a set  of  expansion  ft*nctlons  {M  , n-1,2, . . . ,H},  and  let 


4<W1 


M - T V M 

^ r\Am 


(5) 


where  the  coefficients  V are  to  be  determined.  Substitute  (5)  into  (4) 

n 

and  usf  the  linearity  of  the  H operators  to  obtain 


y V Hf(M  ) + y V H^(M  ) - - 

^ n^  **n  ^ n*<  **n  •«*<: 

n n 


(6) 


Next,  define  a synmetrlc  product 


<A,B> 


B ds 


(7) 


apert. 


and  a set  of  testing  functions  {W  , n-1,2, . . . ,N} , which  may  or  may  not  be 
equal  to  the  expansion  functions . We  take  the  symmetric  product  of  (6) 
with  each  testing  function  W , and  use  the  linearity  of  the  symmetric 
product  to  obtain  the  set  of  equations 

‘'nmtn  “nmtn  mt 

n n 

!if-1,2,.  . . ,N.  Solution  of  this  set  of  linear  equations  determines  the 

coefficients  V and  the  magnetic  current  M according  to  (5).  Once  M 
n -.AV- 

Is known,  the  fields  and  field-related  parameters  may  be  computed  by 

standard  methods. 


[4]  R.  F.  Harrington,  "Field  Computation  by  Moment  Methods,"  Macmillan 
Company,  New  York,  196d. 
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The  above  solution  can  be  put  Into  matrix  notation  as  follows: 

Define  an  admittance  matrix  for  region  a as 

and  an  admittance  matrix  for  region  b as 

iy")  - (.0) 

The  minus  signs  are  placed  In  (9)  and  (10)  on  the  basis  of  power  con- 
siderations. Define  a source  vector 

(11) 

and  a coefficient  vector 

^ * <''Jk«1 

Now  the  matrix  equation  equivalent  to  equations  (8)  is 

[Y^  + (13) 

This  can  be  Interpreted  In  terms  of  generalized  networks  as  two  networks 
[Y®]  and  [Y^]  In  parallel  with  the  current  source  as  shown  In  Fig.  3. 
The  resultant  voltage  vector 

V = [Y®  + (14) 

Is  then  the  vector  of  coefficients  which  give  M according  to  (5). 

It  Is  Important  to  note  that  computation  of  [Y  ] Involves  only  region 
a,  and  computation  of  [Y^]  involves  only  region  b.  Hence,  we  have  divided 
the  problem  Into  two  parts,  each  of  which  may  be  formulated  Independently. 
Once  [Y]  is  computed  for  one  region,  it  may  be  combined  with  [Y]  for  any 
other  region,  making  it  useful  for  a wide  range  of  problems.  For  example, 
the  same  aperture  admittance  matrix  for  radiation  into  half-tnace  would  be 
useful  for  plane-wave  excitation  of  the  aperture,  waveguide  excitation, 
and  cavity  excitation. 


III.  LINEAR  MEASUREHENT 


A linear  measurement  Is  defined  as  a number  which  depends  linearly 

on  the  source.  Examples  of  linear  measurements  are  con^onents  of  the 

field  at  a point,  voltage  along  a given  contour,  and  current  crossing  a 

given  surface.  Measurements  made  In  region  b will  depend  linearly  only  on 

the  equivalent  current  -M.  Measurements  made  In  region  a will  depend 

^ 11 

linearly  on  the  Impressed  sources  ^ ^ , as  well  as  on  the  equivalent  cur- 
rent M.  We  now  Illustrate  these  concepts  with  a particular  example. 

Consider  the  measurement  (computation)  of  a component  H of  magnetic 

m 

field  at  a point  r^  In  region  b.  It  Is  known  that  this  component  can  be 
obtained  by  placing  a magnetic  dipole  K1  ac  r , and  applying  the  reciprocity 
theorem  to  Its  field  and  to  the  original  field  [2,  Sec.  3-8].  The  original 
field  In  region  b Is  given  by  the  solution  to  Fig.  2b.  The  problem  Involving 
the  magnetic  dipole,  called  the  adjoint  problem.  Is  shown  In  Fig.  4.  Appli- 
cation of  the  reciprocity  theorem  to  these  two  cases  yields 


H K£  - - 
mm  j 


M • if'ds 


05) 


apert. 


Here  Is  the  magnetic  field  from  In  the  presence  of  a complete  con- 
ductor, and  H Is  the  component  In  the  direction  of  K£  of  the  magnetic 

n 

field  at  r due  to  -M  In  the  presence  of  a complete  conductor.  To  evaluate 

.vkm  A*- 


(15) , substitute  for  M from  (5)  and  obtain 

H Kl  = > V <-M  ,rf“> 

mm  n n ' 


n 


n n 


(16) 


This  can  be  written  in  matrix  form  as 


H K£  » I“V 
m m 


where  Is  the  transpose  of  a measureraent  vector 


n' 


‘Nxl 


(17) 


(18) 


Note  that  the  elements  of  l!°  are  similar  In  form  to  those  of  given 
by  (11) , except  that  -M  replaces  W . The  minus  sign  difference  reflects 


the  fact  that  the  equivalent  source  in  region  b Is  In  contrast  to 
that  in  region  a which  is  -HM.  Now  substitute  (14)  into  (17)  to  ob- 
tain 

H K1  - -J-  (19) 

in  D 

If  the  magnetic  dipole  is  of  unit  moment,  then  (19)  gives  H at  r 

TW  iiaiH 

directly. 

Every  linear  measurement  in  region  b will  be  of  the  form  (19) . 

For  example,  if  a component  of  E at  r were  desired,  we  would  place  an 
electric  dipole  at  ^ and  apply  reciprocity.  In  general,  a linear 
measurement  involves  applying  reciprocity  to  the  original  problem  and 
to  an  adjoint  problem.  A determination  of  the  sources  of  the  adjoint 
problem  is  a part  of  the  formulation  of  the  problem. 

If  a linear  measurement  is  made  in  region  a,  it  will  involve  a 
contribution  from  the  impressed  sources  added  to  that  from  the 

equivalent  sources  M.  For  example,  instead  of  (19)  we  would  have 

H Ki  - nSa  + f"[Y®  + (29) 

m m m m 

where  is  the  magnetic  field  from  in  the  presence  of  a complete 

conductor.  Also,  in  region  a we  would  define  the  measurement  vector  to 
be 


(21) 


instead  of  (18) , because  the  equivalent  sources  are  4M  in  region  a in 
contrast  to  -II  in  region  b.  Note  that  it  is  the  difference  field  H-^ 
in  region  a (due  to  that  is  directly  analogous  to  the  transmitted 
field  H in  region  b (due  to  -M). 


IV.  TRANSMITTED  POWER 

A quadratic  measurement  is  one  which  depends  quadratlcally  on  the 
sources.  Examples  of  quadratic  measurements  are  components  of  the  Foyntlng 
vector  at  a point,  power  crossing  a given  surface,  cmd  energy  within  a 
given  region.  A particular  quadratic  measurement  of  considerable  interest 
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Is  the  power  transmitted  through  the  aperture,  which  we  now  consider. 


The  complex  power  transmitted  through  the  aperture  is  basically 

- 1 1 E X H**  n ds  (22) 

t I I ^ 

apart. 

where  the  asterisk  denotes  complex  conjugate.  Substituting  from  (1),  we 
have 


P^  - [[  M • H*ds 

t J J 


(23) 


apart. 

This  Involves  only  the  tangential  component  of  ^ which  in  region  b we 
b 

denoted  by  ^(-M).  For  ^we  use  the  linear  combination  (5)  and  obtain 


r (-M)  *=  - y V ^(M  ) 


(24) 


n 


Substituting  this  for  K and  (5)  for  into  (23) , we  obtain 

- y y V V*  [[  M • _^*(M  )ds  (25) 

m n ' ' 


apart. 


If  M are  real,  the  conjugate  operations  can  be  taken  outside  the  inte- 

wQ 

grals.  Moreover,  if  M ■ W (Galerkin's  method),  then  the  negative  of 

the  integrals  in  (25)  are  y“  as  defined  by  (10) , ano' 

mn 


p = y y V V*  y' 

t ^ ^ m n 1 


m n 


* „b* 

nn 


(26) 


This  can  be  written  in  matrix  form  as 


Pj.  *■  v[y'’]*^* 


(27) 


Note  that  this  is  the  usual  formula  for  power  into  network  [Y  ] of  Fig.  3. 


V.  APERTURES  IN  PLANE  CONDUCTORS 

Consider  a conducting  plane  covering  the  z»0  plane  except  for  an 
aperture,  as  shown  in  Fig.  5.  The  two  regions  z>0  and  z<0  are  identical 
half  spaces,  and  hence  their  admittance  matrices  are  the  same.  Therefore, 
we  let 
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Fig.  5.  Aperture  In  a plane  conductor. 


[Y®  + Y^]  - 2[Y^®] 


(28) 

where  [Y*"*]  denotes  the  aperture  admittance  for  the  aperture  opening  Into 
half  space,  say  z>0.  When  the  aperture  Is  covered  by  a conductor,  the 
z«0  plane  Is  a complete  conducting  pleuie,  and  image  theory  appllet-''.  The 
magnetic  current  expansion  functions  are  on  the  surface  of  the  z«0  plane. 
Their  Images  are  equal  to  them  and  are  also  on  the  z>0  plane  [2,  Sec.  3-6]. 
The  result  is  that  [Y*  ] Is  the  admittance  matrix  obtained  using  expansion 
functions  radiating  Into  free  space  everywhere.  This  problem  Is  dual 
to  that  for  the  impedance  matrix  of  a plane  conductor,  a problem  considered 
recently  In  the  literature  [51. 


The  original  excitation  of  the  aperture  Is  by  the  Impressed  sources 
J^,  In  the  region  z<0.  The  In^ressed  field  used  in  the  operator 
equation  (4)  Is  the  tangential  magnetic  field  due  to  , M with  the  aper- 
ture covered  by  a conductor  (Fig.  2a).  In  this  case  the  z“0  plane  Is  a 
complete  conductor,  and  Image  theory  again  applies.  The  result  Is  that 
the  tangential  component  of  jS  over  the  z<°0  plane  when  It  Is  covered  by  a 
conductor  Is  just  twice  what  it  is  for  the  same  sources  In  free  space. 
Hence, 


<*< 


2H*° 


(29) 


where  la  the  tangential  component  of  the  magnetic  field  over  the 

11 

aperture  due  to  the  sources  , M in  free  space.  The  components  of 
the  excitation  vector  defined  by  (11)  are  now 


1^-2  ff  W • J^°ds  (30) 

apert. 

where  H Is  the  mth  testing  function. 

■**Tn 

A case  of  special  Interest  Is  that  of  plane  wave  excitation.  A 
unit  plane  wave  Is  given  by 


[5]  N.  N.  Wang,  J.  H.  Richmond,  M.  C.  Gllreath,  "Sinusoidal  Reaction 

Formulation  for  Radiation  and  Scattering  from  Conducting  Surfaces," 
IEEE  Trans,  on  Antennas  and  Propagation,  vol.  AP-23,  No.  3,  pp.  376- 
382,  May  1975. 
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* r 


(37/ 


ifhere  is  a unit  vector  specifying  the  direction  of  ^ is  the 

propagation  vector  of  magnitude  2ir/A  and  pointing  in  the  direction  of 
propagation,  and  r^ is  the  radius  vector  to  an  arbitrary  field  point. 
These  vectors  are  shown  in  Fig.  5.  The  conq>onents  (30)  of  the  plane- 
wave  excitation  vector  are  then 


m 


W 


(32) 


apert. 

% 

The  synA>ol  has  been  used  for  this  particular  vector  to  distinguish 
it  from  the  more  general  excitation  vector  (30) . 


Similar  simplifications  apply  to  the  adjoint  (measurement) 
problem.  For  the  evaluation  of  a component  of  magnetic  field  at  a 
point  r , a magnetic  dipole  Ki  is  placed  at  the  measurement  point  r . 
This  radiates  in  the  presence  of  a complete  conductor  over  the  z“0 
plane,  and  hence,  analogous  to  (29),  we  have 


ni 


•~t 


(33) 


Here  denotes  the  teingentlal  component  of  H over  the  aperture  from 
K£  when  the  z“0  plane  is  covered  by  a conductor,  and  denotes  that 
from  K£  when  it  radiates  into  free  space.  The  components  of  the  measure- 
ment  vector  defined  by  (18)  are  now 


..m 


JJ 

apert. 


M 


(34) 


where  M is  the  nth  expansion  function. 

-*n 

A case  of  special  Interest  is  that  of  far-fleld  measurement. 

This  is  obtained  by  a procedure  dual  to  that  used  for  radiation  and  scat- 
tering from  conducting  wires  [6].  To  obtain  a component  of  H on  the 
radiation  sphere,  we  take  a source  perpendicular  to  ^ and  let  r^  •+•  ». 


[6]  R.  F.  Harrington,  "Matrix  Methods  for  Field  Problems,"  Proc.  IEEE, 
vol.  55  No,  2,  pp.  136-149,  February  1967. 
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At  the  sane  time  we  adjust  so  that  It  produces  a unit  plane  wave 
In  the  vicinity  of  the  origin.  The  required  dipole  moment  Is  given 


(35) 


and  the  plane-wave  field  It  produces  In  the  vicinity  of  the  origin  Is 


(36) 


Here  u Is  a unit  vector  Ir  the  direction  of  k Is  the  propagation 

vector,  and  r is  the  radius  vector  to  an  arbitrary  field  point.  Again 
these  vectors  are  shown  In  Fig.  5.  The  components  (34)  of  the  far-fleld 
measurement  vector  are  then 


apert. 


u 


(37) 


The  symbol  Is  used  for  this  particular  measurement  vector  to  dis- 
tinguish It  from  the  more  general  measurement  vector  (34) . The  far- 
zone  magnetic  field  Is  now  given  by  (19)  with  given  by  (35) , 

and  the  aperture  admittance  given  by  (28) . Hence 


m Sirr 


(38) 


The  usual  two  radiation  components  H.  and  H.  are  obtained  by  orienting 
Kl  in  the  6 and  if  directions,  respectively. 


A parameter  sometimes  used  to  express  the  transmission  character- 
istics of  an  aperture  Is  the  transmission  cross  section  t.  It  Is  de- 
fined as  that  area  for  which  the  incident  wave  contains  sufficient  power 
to  produce  the  radiation  field  H^  by  omnidirectional  radiation  over  half 
space.  For  unit  Incident  magnetic  field,  this  Is 


T * 


(39) 
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Substituting  froa  (38),  ve  obtain 


T 


32ir 


2 


(40) 


Note  that  t depends  upon  the  polarization  and  direction  of  the  incident 

-i-i 

wave  (via  P ),  and  upon  the  polarization  measured  and  direction  to  the 
measuremeut  point  (via  . 

Another  parameter  used  to  express  the  transmission  characteristics 
of  an  aperture  is  the  transmission  coefficient  T,  defined  as 


T 


trans. 


Inc. 


(41) 


where  P is  the  time-average  power  transmitted  by  the  aperture, 

uirdAS  a 

and  l8  the  free  space  power  incident  on  the  aperture.  The  inci- 

dent power  is 


Inc. 


n S cos  6 


inc 


(42) 


where  n ■ /u/e  is  the  intrinsic  Impedance  of  free  space,  S is  the 
aperture  area,  and  is  the  angle  between  ^ and  n.  The  trans- 

mitted power  is 


P 

trans. 


Re(Pj.) 


where  Re(P^)  denotes  the  real  part  of  P^,  given  by  (27)  with 
[Y^]  = [Y^®].  Hence 


T 


1 


n S cos  9 


Inc. 


Re(^[Y^®]*^*) 


(43) 


(44) 


.’lote  that  T depends  on  both  the  direction  of  incidence  and  oii  the 
polarization  of  the  incident  wave. 

Finally,  because  of  symmetry  about  the  z*0  plane,  the  difference 
field  which  exists  in  the  region  z<0  is  slnq;>ly  related  to  the 
transmitted  field  which  exists  in  the  regnn  z>0.  The  difference  field 
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in  the  region  z<Q  is  produced  by  an  equivalent  current  M on  a plane 
conductor  over  the  z-0  plane.  By  image  theory,  it  is  also  the  field 
produced  in  the  region  z<Q  by  the  source  2M  in  free  space.  As  noted 
earlier,  the  transmitted  field  in  the  region  z>0  is  produced  by  the 
source  -2M  in  free  space.  Hence,  the  difference  field  in  the  region 
z<0  and  the  negative  of  the  transmitted  field  in  the  region  z>0  are 
botli  produced  by  the  same  magnetic  current  ^ radiating  In  free  space. 
In  other  words,  the  difference  field  is  the  negftive  of  the  analytic 
continuation  of  the  transmitted  field  into  the  region  z<0.  This 
applies  to  electric  fields  as  well  as  to  magnetic  fields. 


VI.  WAVEGUIDE-FED  APERTURES 

Consider  now  a uniform  waveguide  feeding  an  aperture  in  a con- 
ducting plane,  as  shown  in  Fig.  6.  In  general,  the  aperture  may  be  of 
different  size  and  shape  than  the  waveguide  cross  section.  The  half- 
space region  z>0  is  the  same  as  in  the  previous  problem.  Fig.  5,  and 
the  analysis  of  the  preceding  section  applies.  An  analysis  of  the 
waveguide  region  Is  given  below. 

Let  the  excitation  of  the  waveguide  be  a source  which  produces  a 
single  mode,  of  unit  amplitude,  incident  on  the  aperture.  This  mode 
(usually  the  dominant  mode)  is  denoted  by  the  index  o.  The  field  tan- 
gential to  the  z-dlrectlon  can  then  be  expressed  In  modal  form  as 
[2,  Sec.  8—1] 


E 

-t 


-y  z 

f.  ° e + y r.e 
-o  J 1 


Yi^ 


St 


-y  z 
V O 
Ye  u 

o ~z 


e - y r.Y,  e 

-'O  “11 

1 


YiZ 


u 


(45) 


It  Is  assumed  that  all  modes,  TE  and  TM,  are  Included  in  the  summation. 
The  sre  modal  propagation  constants 


- jk 

(XA^'^ 

X < X 

(X^/X)^ 

X > X 

(46) 
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Fig.  6.  Waveguide-fed  aperture  In  a conducting 
plane. 


where  is  the  1-th  mode  cut-off  wavelength,  and  ■=  2u/Xj^  is  the 
i-th  code  cut-off  wavenumber.  The  are  the  modal  characteristic 
admittances 


V 

"1 


TE  modes 


TM  modes 


(47) 


F Is  the  reflection  coefficient  for  the  o-th  mode,  and  F,  Is  the 
o 1 

complex  amplitude  of  the  -z  traveling  component  of  the  1-th  mode. 

The  e^  are  normalized  modal  vectors,  so  that  the  modal  orthogonality 
relationships  are 


■r 

. 

guide 


Si 


i * j 
1 » j 


(48) 


where  the  Integration  is  over  the  waveguide  cross  section. 

To  evaluate  the  aperture  admittance  (9)  In  the  waveguide  region, 
we  consider  a single  expansion  function  ^ on  the  z*0  plane  In  the 
waveguide  region.  The  tangential  field  produced  by  will  be  of  the 
form  (45),  except  that  there  Is  no  Incident  wave.  Hence,  this  field 
Is 


■ p„l  ' 

■ - 1 a. " a 


(49) 


where  the  are  modal  amplitudes.  At  z*0  we  have 


M “ u X E. 
Ml 


z-0  1 ^ 


(50) 


Multiply  each  side  of  this  equation  scalarly  by  u^  x and  Integrate 
over  the  waveguide  cross  section,  obtaining 

T f f 

M • u X e.ds  "•yA.  (u  xe)«(u  x e.)ds 

"nljj  -^*1 


guide 


guide 


(51) 


By  orthogonality  (48),  all  terns  of  the  summation  are  zero  except  the 
l*j  tern.  Hence, 


If 


M 

•n 


u X e . ds 
HL 


(52) 


apert. 

We  have  replaced  the  Integral  over  the  waveguide  cross  section  by  one 
over  the  aperture,  since  M exists  only  In  the  aperture  region.  The 
elements  of  the  aperture  admittance  matrix  (9)  are  nov  given  b> 


r 


& 


yWg 

sn 


ff  W • ^(M  )ds 
apert. 


(53) 


where  the  superscript  wg  denotes  waveguide.  The  H of  (53)  Is  given  by 
the  second  equation  of  (49)  evaluated  at  z>0,  or 


jWg 

m 


iw.  If 


w 

MO 


u X e.ds 


apert . 


Now  define  the  constants 


ml 


If 


W 

<%»yp 


u X e , ds 
~z 


apert. 


(54) 


(55) 


which  are  similar  In  form  to  the  of  (52) . The  elements  (54)  then 
are  given  by 


mn 


A .Y, 
nl  ml  1 


(56) 


Hence,  all  elements  of  the  waveguide  admittance  matrix  [Y^'^]  are  linear 

combinations  of  the  modal  characteristic  admittances  Y^.  For  Galerkln's 

method,  W > M and  the  A and  B . are  equal. 

-.n  ~n..  *>nl 

We  next  evaluate  the  equ  valent  magnetic  current  M,  given  by  (5). 
The  Incident  field  Is  given  by  the  first  term  on  the  right-hand  side 
of  (45).  When  the  aperture  Is  covered  by  a conductor,  the  waveguide 
Is  terminated  by  a conducting  plane.  According  to  Image  theory,  the 
tangential  magnetic  field  at  z»0  Is  then  just  twice  the  Incident  wave, 
or 


) 
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H"  - 2f  u X e (57) 

^ 0'**Z  **0 

This  is  the  used  in  (11)  to  evaluate  the  excitation  vector  Hence, 
the  components  of  the  excitation  vector  are 


- 2Y  f I W ‘ u X e ds 
m o J J '■*®  •** 


2Y  B 

O QO 


apert. 

The  total  aperture  admittance  matrix  is 

[Y®  + y'’]  - [y''*  + Y^®]  (59) 

where  [y'^®]  is  the  waveguide  aperture  admittance  and  [Y^®]  is  the  half- 
space aperture  admittance.  The  coefficient  matrix  V is  given  by  (14) 
with  the  admittance  matrix  given  by  (59) , or 

i)'  > (60) 

Finally,  the  equivalent  magnetic  current  M is  given  by  (5)  where  the 
coefficients  are  the  components  of  V, 

Once  M is  found,  the  modal  amplitudes  F.  in  (45)  can  be  evaluated 
from  (1)  and  the  orthogonality  properties  of  the  modes.  From  (1)  and 
(45) , we  have 


M=u  xE^  “U  xe  + Tr,u  xe  (61) 

— «-z  »o  “ l~z  ^ 

z”0  i 

Now  multiply  each  side  scalarly  by  u ^ ^ and  Integrate  over  the  wave- 
guide  cross  section.  By  the  orthogonality  relationships  (48) , all  terms 
of  the  summation  vanish  except  the  term  l->j.  The  result  is 


• f 

M • u 


u X e,ds 
~z  1 


guide 


] + r i - 0 

o 


1 y 0 


Here  the  integration  over  the  guide  can  be  changed  to  that  over  the 
aperture  because  M * 0 except  in  the  aperture.  Substituting  for  II 
from  (5)  into  (62),  and  using  the  de:^initlon8  (52),  we  have 


X 


(63) 


V A - 1 + r 

n no  o 


yvA=r, 

^ n nl  i 
n 


Finally,  by  defining  nodal  measurement  vectors  as 


(64) 


and  using  (60),  we  can  write  (63)  as 


1 + r - A Iy''®  + (65) 

o o 

and,  for  1 0, 

Ff  - A^[y''®  + (66) 


The  parameter  of  most  Interest  Is  T , the  reflection  coefficient  of  the 

o 

Incident  mode.  This  Is  often  expressed  In  terms  of  an  admittance 


1 - r 


1 + r o 

o 


(67) 


which  Is  the  equivalent  aperture  admittance  seen  by  the  Incident  mode. 

The  region  ^>0  for  waveguide-fed  apertures  Is  the  same  half-space 
region  as  existed  in  the  previous  problem  of  an  aperture  In  a conducting 
plane.  Hence,  evaluation  of  the  fields  in  terms  cf  M In  this  region  Is 
done  the  same  way  as  in  Section  V.  For  example,  the  u cooiponent  of  the 
far-zone  magnetic  field  at  a point  r is  given  by 


H 

m 


^ p"[y"«  + Y>"1-¥ 

m 


(68) 


which  is  (38)  with  the  term  [2Y^®]  replaced  by  (60).  The  excita- 
tion vector  has  elements  given  by  (58),  and  the  far-fleld  measurement 
vector  ^ has  elements  given  by  (37).  The  power  gain  pattern  is  the  ratio 
of  the  radiation  Intensity  in  a given  direction  to  the  radiation  intensity 
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which  would  exist  if  the  total  power  were  radiated  uniformly  over  half 
space,  or 


G 


Re(P^) 


(69) 


Subp«-ituting  for  H from  (68),  we  have 
m 


^ ~ 8rke(P  ) (70) 

where  P^  is  given  by  (27).  Note  that  this  gain  is  a function  of  the 
H component  measured,  as  well  as  direction  to  the  field  point. 


VII.  DISCUSSION 


A genera]  iormulaticn  for  aperture  problems  using  the  method  of 
moments  has  b»en  given.  The  solution  is  expressed  in  terms  of  aperture 
admittance  matrices,  one  fo.:  each  region.  While  the  exposition  has 
been  written  in  terms  of  a single  aperture  the  formulation  applies 
equally  well  to  multiple  apertures.  The  extension  of  the  theory  to 
several  regions  is  also  relatively  straightforward.  An  example  of 
this  more  general  case  ic  many  waveguides  radiating  into  half  space,  a 
problem  treated  bjf  Galindo  and  Wu  [3].  Another  example  is  a coaxial  line 
feeding  a cavity,  which  in  turn  radiates  into  half  space  through  an  aper- 
ture. In  this  case  the  generalized  network  equivalent  would  be  a network 
representing  the  cavity,  connected  through  some  ports  to  a network  repre- 
senting the  coaxial  line,  and  connected  through  other  ports  to  a network 
representing  half  space. 

Explicit  formulations  have  been  given  for  two  problems,  that  of  an 
aperture  in  a conducting  plane  with  plane-wave  excitation,  and  that  of  a 
waveguide  feeding  an  aperture  in  a conducting  plane.  General  computer 
programs  have  been  developed  for  rectangular  apertures  in  conducting 
planes,  and  for  rectangular  apertures  in  rectangular  waveguides.  These 
programs  will  be  described  in  future  reports.  It  is  also  planned  to 
apply  the  general  theory  to  develop  computer  programs  for  rectangular 
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apertures  In  rectangular  cavities.  We  will  then  have  programs  which  can 
be  used  for  computing  the  behavior  of  rectangular  apertures  coupling  two 
regions  which  are  any  permutation  of  these  three  cases.  Other  applica- 
tions of  the  general  theory,  such  as  to  apertures  and  cavities  of  arbitrary 

I 

shape,  and  to  multiple  aperture  problems,  are  contemplated  for  future  work. 
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